The Wardle Instability in Interstellar Shocks: I. Nonlinear Dynamical Evolution

The nonlinear evolution of unstable C-type shocks in weakly ionized plasmas is studied by means of time-dependent magnetohydrodynamical simulations. This study is limited to shocks in magnetically dominated plasmas (in which the Alfvén speed in the neutrals greatly exceeds the sound speed), and microphysical processes such as ionization and recombination are not followed. Both two-dimensional simulations of planar perpendicular and oblique C-type shocks, and fully three-dimensional simulation of a perpendicular shock are presented. For the cases studied here, the instability results in the formation of dense sheets of gas elongated in the direction of shock propagation, and oriented perpendicular to the magnetic field. The formation of a weak J-type front is associated with the growth of the instability from an equilibrium shock structure. After saturation, the magnetic field structure consists of arches which bow outwards in the direction of shock propagation, and which are anchored by the enhanced ion-neutral drag in the dense sheets. In analogy to the magnetic buoyancy (Parker) instability, saturation occurs when the magnetic tension in the distorted field lines is balanced by drag in the sheets. For the magnetically dominated shocks studied here, the distortions in the magnetic field which produce saturation are very small. Nonetheless, the enhancements of the ion and neutral densities in the sheets is very large, between two and three orders of magnitude compared to the preshock values. At these high densities, recombination processes may be important. The sheets evolve slowly in time, so that shocks propagating in a homogeneous medium may leave behind a network of intersecting filaments and sheets of dense gas elongated in the direction of shock propagation and perpendicular to the mean field. The temperature structure and emission properties of unstable C-type shocks in the nonlinear regime are presented in a companion paper (Neufeld & Stone 1997). Subject headings: shock waves instabilities MHD methods: numerical ISM: kinematics and dynamics, magnetic fields